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0^ ' The momentum dependence of the mean-field contribution to the real part of the optical model 

potential is investigated employing realistic nucleon-nucleon interactions. Within a non-relativistic 
' approach a momentum dependence originates from the non-locality of the Fock exchange term. De- 

ducing the real part of the optical model from a relativistic Dirac Brueckner Hartree Fock approx- 
imation for the self-energy of the nucleons yields an additional momentum dependence originating 
from the non-relativistic reduction of the self-energy. It is demonstrated that large Fock terms are 
required in the non-relativistic approach to simulate these relativistic features. A comparison is 
made between a local density approximation for the optical model and a direct evaluation in finite 
nuclei. 
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CN ■ I. INTRODUCTION 

o ■ 

QS\ ■ The nucleon - nucleus optical potential is a very important tool for the analysis of nuclear reaction data. Therefore 
a lot of empirical information on nucleon - nucleus scattering for a large spectrum of energies and various nuclei has 
been collected and the global fits to these optical potentials provide a very condensed information on nuclear 
Q ■ structure properties. It is obvious that theoreticians have devoted much interest to the study of this important 
^ . quantity. Rather than trying to give a short overview on such investigations we refer to some recent review articles 
on this subject 

From a theoretical point of view the real part of the optical potential can be split into a mean-field contribution 
' and a dispersive correction ||^ . The name dispersive contribution originates from the fact that this real part is related 
to the imaginary part of the optical potential by a dispersion relation. This dispersion relation yields an explicit 
energy-dependence also for the real part. In principle, the dispersive contribution could be determined from the 
empirical data on the imaginary component. This would require the knowledge of the imaginary part at all energies. 
For large energies, however, the imaginary part gets contributions from nucleonic excitations, due to the strong 
short-range components of a realistic NN interaction and the excitation of sub- nucleonic degrees of freedom, like 
particle production, which is beyond the scope of the present discussion. Therefore applications of such an dispersive 
optical-model analysis typically consider only the imaginary part at low energies ]lOt and use some kind of subtracted 
dispersion relation ||]. 

If one tries to describe the remaining mean-field contribution of the optical potential in terms of a local potential 
to be used in a Schroedinger equation, one observes a significant dependence of this mean-field part on the energy or 
asymptotic momentum of the scattered nucleon There are two possible sources for this momentum dependence. 

Within a non-relativistic approach one may attribute this momentum dependence to the non-locality of the exchange 
part in the Hartree-Fock (HF) or Brueckner-Hartree-Fock (BHF) approximation for the mean-field part of the optical 
potential. Transforming this non-local HF potential to an equivalent local potential leads to such a momentum 
dependence 0,0 (^^^ chapter 2). 
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Another source for the momentum-dependence of the mean-field part of the optical potential is of relativistic 
origin. Relativistic approaches have been very successful to describe bulk properties of nuclei. This is true for 
more phenomcnological models, like the a — uj model of Walccka and Serot [Q, but also for microscopic Dirac-BHF 
calculations based on realistic nucleon-nucleon (NN) forces [14-1^. Characteristic for all these calculations is the 
relativistic structure of the nucleon self-energy, which contains a large Lorentz-scalar and a large time-like vector 
component. In a simple Dirac-Hartree calculation these components are local and energy independent. Reducing 
the two coupled Dirac equations for the upper and lower component to one Schroedinger equation, one derives from 
these scalar and vector components a central Schroedinger equivalent potential (SEP), which depends on the energy 
or momentum of the nucleon |p^ . 

The scalar and vector components of the relativistic self-energy obtained in a Dirac-HF or Dirac-BHF calculation 
are non-local. Therefore both sources for the momentum dependence of the mean-field part of the optical model 
discussed above get relevant. It is one aim of the present investigation to explore the importance of these two effects 
and their interplay. This study is based on the self-consistent Dirac-BHF calculations, which were recently performed 
directly for finite nuclei ||l^,|2^ . 

There have been many very successful attempts to deduce the mean field part of the optical potential directly from 
the NN T-matrix or G-matrix [Q within a non- relativistic [^l|-p4| and relativistic framework |^5|-|27|]. The impulse 
approximation which is used in most of these investigations is rather successful in particular for proton scattering at 
energies above 100 MeV. Medium effects are typically derived from nuclear matter and taken into account in a local 
density approximation [^ . 

The main emphasis of the present investigation is to study the differences and similarities in non-relativistic and 
relativistic approaches in particular at lower energies. These energies are important for the extrapolation of the single- 
particle potential at negative energies, connecting the optical potential to the binding properties. For that purpose 
we consider simple parameterizations of G-matrix approaches, derived from realistic NN forces. In the case of the 
non-relativistic approach this will be the so-called M3Y parameterization . In the case of the relativistic treatment 
this will be the effective meson parameterization of 19 2^]. This enables us to study and compare the origin of the 
momentum dependence obtained in the non-relativistic and relativistic approach in detail. Furthermore we would like 
to compare the results for nucleon-nucleus scattering obtained from a non-local Dirac self-energy which is calculated 
directly for the finite nuclei to results which are obtained in the frequently used local-density approximation. 

After this introduction we will briefiy recall the connection between the non-locality and the momentum dependence 
of an equivalent local mean field in section 2. In section 3 we will shortly describe the momentum dependence of the 
SEP derived from a local Dirac self-energy. There we will also describe an efficient method to calculate scattering 
phase shifts from a non-local relativistic self-energy. The results of our investigation are presented and discussed in 
section 4 and the conclusions are summarized in a final section 5. 



II. NON-LOCALITY AND EXCHANGE EFFECTS 

More than 15 years ago Bertsch et.al. determined a very efficient parameterization of an effective NN interaction 
from the G-matrix evaluated for realistic forces |2^]. These so-called M3Y potentials are given in terms of local 
Yukawa potentials. The central part of the interaction can be written as a function of the distance r between the 
interacting nuclcons 

^M3Y ^ ^o(^) + ^rir)fi ■ T2 + Vair)(Ji ■ (72 + KrOi • 0^2 Tl ■ . (1) 

The operators Ti and Gi denote the Pauli matrices acting on the isospin and spin of nucleon i, respectively. The radial 
shapes of the functions Va {r) are given in terms of 2 Yukawa potentials 

V.ir)^tc^^^-^^^^. (2) 
fe=i ^'"'^ 
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only for the spin-isospin channel Var a third Yukawa term has been added, which represents the central component 
of the pion-exchange potential 

^ ^ ar — ^ . . o ''''71 



3 47r Am? fi^r 
= C^'E^hl^ . (3) 

Here — m-^/{hc) — 0.7072 fm~^ and C-^ = 3.4877 MeV, which corresponds to a pseudoscalar coupling constant 
for the pion of 5^/(47r) w 13.9. Assuming for the mass parameter /i^ in Eq.(||) values of 4 fm~^ and 2.5 fm^^, for 

(k) . 

k = 1 and 2, respectively, the 8 strength parameter Ca were adjusted to reproduce G-matrix elements of the Reid 
potential and the Elliot potential |31| in an oscillator basis. The resulting parameters are listed in table 1. 

For the mean-field contribution of the optical potential for scattering of nucleons on closed shell nuclei with N=Z 
only the spin- and isospin averaged part of the NN interaction is relevant. For the direct interaction leading to the 
Hartree part of the potential this is identical to the scalar isoscalar term in Eq.(|l|) 

Vair = Mr) 

= 7999.0^^^4^ - 2134.25^^4^ ■ (4) 
4r 2.5r 

The resulting direct contribution to the optical potential is given by 

UdiA^) = J d'r'p,ir')V^-^i\r~ r^\) (5) 
with the nuclear density Pb(r) — p{r,r) obtained in the Hartree- Fock approximation 

N 

p(r,r')=^<l'*(r)$.(r'). (6) 



Here $i stands for the single-particle wavefunctions occupied in the Hartree-Fock approximation for the target nucleus. 
The exchange part of the NN interaction is evaluated as the scalar isoscalar part of the Fierz transformation of 
Eq.(|l]) and is given as 

Vexc [V^{r) + Wr{r) + W„{r) + W^r{r)] 

= 4631.375^^P(^ - 1787.125^^Pi^ - 1 MlA^'-^^t^^^ . (7) 
4r 2.5r 0.7072r ^ ' 

This local parameterization of the G-matrix and slight modifications of it have been used for various nuclear studies 
including the evaluation of the optical potential for heavy-ion scattering . The exchange part of the NN interaction 
yields the Fock contribution to the mean field which is non-local 

t/exc(r,r') - p{r,r')Ve-Kc{\r^r\) (8) 

The local equivalent potential is defined by the requirement 

C4xcM*i^M = / d'r'UexcM^E^r') , (9) 

where ^'^(r) is the solution of the Schroedinger equation for a nucleon at energy E moving in the single-particle 
potential U = U^y^^ + J7exc- Negele and Vautherin |Q demonstrated that Eq.(^) can well be approximated by 

C/^xcM = / d'r'Vex^{\r-A)p{ryf^^^^^^^ , (10) 
with the local approximation to evaluate k 
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k'{r)^-^[E-Uir)] . (11) 

h 

It is evident from Eq.(^ that the local equivalent potential depends on the energy E or the momentum k{oo) of 
the scattered nucleon. In Eq.([l0|) this momentum-dependence is hidden in the local momentum fc(r), which must be 
calculated in a self-consistent way from Eq.(|Tl|). 

The mixed density p(r, r') entering the exchange part of the mean field potential is frequently approximated by the 
so-called Slater approximation obtained in nuclear matter 

with the Spherical Bessel function ji (x) | |35[ | and the local Fermi momentum, which is related to the baryon density 
Pb at the central point between r and r' by the nuclear matter relation 

4 = (13) 



III. DIRAC EQUATION AND SCHROEDINGER EQUIVALENT POTENTIAL 

A simple parameterization of Dirac-Brueckner-Hartree-Fock (DBHF) calculations for nuclear matter employing 
realistic NN interactions has been presented in Ref. ||l^,^. For each density pf, two coupling constants for a scalar 
(Gcr) and vector meson (G^^) exchange have been determined such that a simple Dirac Hartree calculation within 
the a — Lu model reproduces the bulk properties of nuclear matter at this density obtained from DBHF. The density 
dependence of the coupling constants reflects the density dependence of the correlations taken into account in DBHF. 

Using these density dependent coupling constants one can determine the nucleon self-energy in a Dirac-Hartree 
(DH) approximation 

^DH = SbH-7°5:5bH (14) 
with a scalar and vector contribution given by 

Sf,rr(r) = -G,(r)m, / r'^dr'G„{r')ps{r')Io{m„r<)Ko{m,ry) (15) 
Jo 

^Dh(^) = -GUr)m^ / r'^dr'G^{r')pb{r')Io{m^r<)Ko{m^ry) (16) 







In these equations the coupling constants G^ (r) and G^j (r) denote the density dependent coupling constants discussed 
above calculated at the nuclear density which is equal to the baryon density pi, of the nucleus at the position r. The 
meson mass parameters were chosen to be TOct=550 MeV and mi^=78S MeV in accordance with the realistic OBE 
potentials fl^ . The functions and Ki^(x) arise from the multipole expansion of the meson propagator in 

coordinate space and are defined using the modified spherical Bessel functions / and K ]3^ ]: 

Il{x) = j= — Kl{x) = j= . (17) 

y/x yjx 

The baryon density and the scalar density ps are defined by 

a 

^^(^) - i E(2J« + 1) [.9a(0 - fair?] , (18) 
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with a summation index a ranging over all occupied shells. The functions ga and fa refer to the small and large 
components of the Dirac spinor 

< r\a vl/„(f) - f ^ /;^''^ „ 1 y.^.,mA^) , (19) 

using the notation of Rcf. . The Dirac spinors are obtained by solving the Dirac equation including the nucleon 
self-energy defined in Eq. (|l^) . The definition of the self-energy and the solution of the Dirac equation are connected 
by a self-consistent requirement. Note, however, that in the present model, the baryon density pb obtained from 
solving the Dirac equation enters the definition of the self-energy not only as displayed explicitly in Eqs.(^5|) and jl^ ) 
but also via the density dependent coupling constants. 

The contributions to the self-energy S^jj and Sj-^jj are local in the Dirac Hartree approximation. The Dirac 
equation corresponds to two coupled equations for the large (g) and small (/) component of the resulting Dirac spinor. 
One can eliminate the small component in a standard way and obtain a Schroedinger equation for a wavefunction 
which at large distances agrees with the large component of the Dirac spinor [|8|j3^ . The central part of the potential 
in this Schroedinger equation, called the central part of the Schroedinger Equivalent Potential (SEP) can be written 

TO 2m 



with 

2 



3 

'-^Darwin (^) ~ ^ 



1 dD{r) 
D{r) dr 



1 dD{r) 1 d^D{r) 



rD{r) dr 2D{r) dr"^ 



D{r) = m + E + ^'{r) + Y.°{r) . (21) 

Therefore we see that this SEP depends on the energy of the scattered nucleon E = — m or its asymptotic 
momentum. This is true already if we employ the Dirac - Hartree approximation, i.e. deal with local functions for 
the scalar and vector component of the nucleon self-energy. 

If the self-energy is evaluated in the Dirac-Hartree-Fock approximation, the contributions are non-local. In order to 
evaluate the scattering observables obtained from a solution of the Dirac equation with such a non-local self-energy, we 
use the technique already described in Ref. ||2^. In this method the radial functions ga{r) and fair), which define the 
stationary solutions of the Dirac equation (see Eq.([l9|)) are expanded in a discrete basis of spherical Bessel functions. 
The wave numbers for this basis are chosen such that this discrete basis is a complete orthonormal basis in a sphere 
of radius D. This radius is chosen to be large enough that the results for the wavefunctions obtained at positive 
energies show the asymptotic form at r ^ D. A typical value for 13 is _D = 20 fm for the light nuclei being considered 
here. With this expansion for the wavefunctions the Dirac equation is rewritten in form of an eigenvalue problem and 
the eigenvalues (Ea) and eigenvectors (the expansion coefficients for ga and fa) are determined by a simple matrix 
diagonalisation . The expressions for the matrixelements of this matrix to be diagonalized are published elsewhere 

Due to the choice of the basis for the expansion (Bessel functions, which vanish at r = D: ji{knD) = 0) also the 
wavefunctions resulting from this diagonalisation must vanish at r = D. This implies that we obtain results for the 
bound single-particle states and for those selected energies of scattering states E,^ for which the wavefunction vanishes 
at r = D. Assuming that the wavefunction has reached its asymptotic form, see discussion above, we can deduce the 
phase shifts 6i at these selected energies from the relation 

g,iD) = cos5iji{KD)^sinSiyi{KD) = 0, (22) 

where yi stands for the Neumann function and the asymptotic wavenumbcr fc,^ is related to the energy E^ by 



E^ ^ E^ + m ^ ^Jkl + to2 . (23) 
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This method for the calculation of scattering phase shifts is a relativistic extension of the scheme employed in Refs. 
p8| , ^ for non-local potentials in a Schroedinger equation. 

Also in the case of the Dirac-Hartree-Fock approximation one can try to find a local potential which used in a 
Schroedinger equation yields the same result for the scattering phase shifts as obtained from the Dirac-Hartree-Fock 
approximation. It is evident that this SEP will exhibit an energy dependence due to the Dirac effects but also due to 
the non-local Fock exchange terms as discussed in the previous section. 



IV. RESULTS AND DISCUSSION 



A parameterization of the Dirac-Brueckner-Hartree-Fock G-matrix calculated in nuclear matter for the version A 
of the Bonn OBE potential has been presented in Ref. |2^]. Employing this parameterization of the G-matrix 
in terms of a density-dependent a and w meson exchange in Dirac-Hartree calculation of finite nuclei, results were 
obtained for the binding energy and radius of closed shell nuclei like ^^O and ^"Ca, which were in fair agreement 



with the experimental data |19 40[|. The baryon densities and scalar densities obtained in such calculations for ^^O 
and '^^Ca are displayed in figures 1 and 2, respectively. These density distributions will be the starting point of our 
further investigations. 

As a first approach we consider the scalar Sj^jj and vector S^jj contribution to the nucleon self-energy in the 
Dirac-Hartree approximation, which are related to these densities via the relations Eqs.(p^) and (p^. From these 
local scalar and vector contributions one can calculate a local SEP as discussed in the preceding section (see Eq. (^o|) ) . 
These local SEP are displayed in the lower parts of figures 1 and 2 for two different energies of the scattered nucleon. 
The attraction of the SEP for nucleons with an energy of 100 MeV is reduced by more than a factor 2 as compared 
to the SEP for nucleons at an energy of 5 MeV. 

This energy dependence of the SEP, which is of course well known, is very important to obtain predictions for the 
elastic scattering which are in agreement with the experimental data. As an example we show in the upper parts of 
figures 3 and 4 results for the I = phase shifts for elastic scattering of nucleons at various energies on ^^O and ^'~'Ca. 
As a bench-mark, representing the experimental data, we also show in these figures the phase-shifts obtained from 
a fit to the experimental data Q. An acceptable agreement with these "experimental" values for the phase shifts 
is obtained, if we evaluate the phase shifts in the Dirac-Hartree approximation or, equivalently, consider the SEP 
for each energy consistently. If, however, the SEP calculated for a nucleon energy of 5 MeV would be used for the 
scattering of nucleons at higher energies as well, one obtains a substantial deviation (see dashed curves in the upper 
parts of figures 3 and 4) . 

As an alternative approach we now consider the non-relativistic parameterization of the G-matrix provided by the 
M3Y approach [ p9[ . It should be mentioned, however, that other non-relativistic parameterizations | ]32[ | show similar 
features. As we want to explore the differences between a relativistic and non-relativistic evaluation of the mean- 
field contribution to the optical potential, independent on the underlying baryon density, we use the same density 
distribution as in the Dirac-Hartree approach. For these densities we determine the scattering potential, including the 
exchange term in the local form of Eq. ( p^ with the Slater approximation of Eq.(p^) for the mixed density p{r, r'). 

The resulting local potentials are also displayed in figures 1 and 2 (dashed curves) again for nucleons with energies 
of 5 and 100 MeV. The shapes of the M3Y potentials deviate from those obtained in the Dirac-Hartree approach 
in a distinct manner. The M3Y potentials tend to be less attractive at the surface and more attractive close to the 
center. The energy dependence of these local potentials, however, seems to be quite similar to the one obtained in 
the Dirac-Hartree approach. 

This is supported by the inspection of the scattering phase shifts in the lower parts of figures 3 and 4. Also for 
the M3Y potentials we obtain an acceptable agreement with the "experimental" data if the energy dependence of the 
local potential is taken into account. If, however, a scattering potential determined for an energy of 5 MeV is used 
for the scattering of nucleons with higher energies as well, the calculated phase shifts are far above the experiment at 
energies around 100 MeV. 
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The results obtained for the phase shifts display a remarkable similarity comparing the Dirac-Hartree and the 
M3Y approach. Similar results are also obtained for other angular momenta and the total cross section. For angular 
momenta I different from zero, and consequently also for the total cross section, the Dirac-Hartree approach tends to 
give slightly larger results. This can be related to the larger radii of the Dirac-Hartree potentials. 

At this stage one may expect that a combination of the energy dependence in the local potential due to the Dirac 
effects and due to the Fock exchange effects, as considered in a Dirac-Hartree-Fock calculation should lead to an 
energy dependence of the calculated phase shifts, which is larger than in the Dirac-Hartree or in the non-relativistic 
Hartree-Fock approximation and therefore exceeds the energy dependence of the empirical data. 

This, however, is not the case. To demonstrate this, we use the Dirac-Hartree-Fock parameterization of the G- matrix 
evaluated in nuclear matter again for version A of the OBEP (including tt and density dependent a and uj exchange). 
We take the same density distributions for ^^O and ^°Ca as discussed above. (The self-consistent Dirac-Hartree-Fock 
densities deviate from the Dirac-Hartree densities used here and are not considered as we would like to study the 
differences in the self-energies which are not due to different densities.) For these ingredients one can calculate the 
non-local Dirac-Hartree-Fock self-energy and evaluate the phase using the technique discussed in section 3. 

Results for the / = phase shifts are displayed in figure 5 for the case of the nucleus ""^Ca. These phase-shifts are 
rather close to those obtained in the Dirac-Hartree approximation. In fact the energy dependence obtained in the 
Dirac-Hartree-Fock approximation is even slightly weaker than the one obtained in the Dirac-Hartree calculation. 

Figure 5 also shows results for the scattering phase shifts calculated in a local density approximation (LDA) 
directly from the nucleon self-energy in nuclear matter ||2^ . This means that one determines the different components 
of the self-energy for various densities of nuclear matter and identifies the local self-energy at position f with the 
corresponding nuclear matter result calculated for the baryon density at this position. For the scalar part of the 
self-energy this means 

SLDA(-'^) = %uc.Mat.(^M'^) (24) 

where Mat ('^' scalar part of the nucleon self-energy calculated in nuclear matter of density p for 

a nucleon with momentum fc, which is equal to the asymptotic momentum of the scattered nucleon. Equivalent 
expressions define the other components of the self-energy. Inserting this expression into Eq. ( pO| ) one obtains the 
SEP in this local density approximation. The phase shifts can be calculated either directly from the solution of the 
Dirac equation with a self-energy according to Eq. (p^ ) or by solving a Schroedinger equation for the SEP. 

From figure 5 we see that the LDA yields results which are in good agreement with the phase shifts calculated 
in Dirac-Hartree and Dirac-Hartree-Fock approximation. Very similar agreement is also obtained in the case of the 
nucleus ^^O, which is not presented here. In particular the energy or momentum dependence of the SEP seems to 
be contained already in the nuclear matter approximation. This can also be seen from the plot of the SEP in figure 
2. The LDA yields predictions for the local equivalent potential, which are in very good agreement with the Dirac- 
Hartree results. In particular the energy dependence is well reproduced. The radial shapes deviate in a characteristic 
way. The LDA reflects the fluctuations in the baryon density more drastically and exhibits a smaller radius than the 
Dirac-Hartree approximation. This is due to the suppression of finite range effects in the nucleon-nucleon interaction, 
which is characteristic for local density approximations . 

From this comparison one can see that the discussion of nuclear matter results using a local density approximation 
yields valuable information on general features of finite nuclei, like the momentum dependence of the optical model. 
It must be kept in mind, however, that the local-density approximation in general is not very reliable in predicting 
specific observables with high accuracy [Q. As an example we would like to mention that the LDA overestimates 
the calculated binding energy per nucleon by 3.2 MeV and 2.3 MeV for ^^O and *°Ca, respectively. Furthermore it 
yields smaller radii for the ground states of these nuclei than the corresponding direct calculation for the finite nuclei 
(see also figure 2 and discussion of the SEP above). 

In our investigation we find that the energy- or momentum-dependence of the optical potential can be described 
within a non-relativistic as well as in a relativistic framework. This is the same situation as in the empirical approaches 
to describe nucleon- nucleus scattering Q . The physical origin of the momentum-dependence of the local equivalent 
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potentials displayed in figures 1 and 2 is quite different in the two approaches. Within the non-relativistic description, 
the momentum dependence is determined by the strong Fock exchange term in the M3Y parameterization. If one 
calculates the binding of nuclear matter with this effective interaction (see Eq.(^, one finds that the dominating 
contribution to the potential energy arises from the Fock exchange terms. This is displayed in figure 6, where the 
direct Hartree contribution to the potential energy of nuclear matter is compared to the total Hartree-Fock result for 
various densities. One finds that the Hartree term only yields about 25 percent of the total attraction. 

These strong exchange terms are required to provide a sufficient momentum dependence for the optical potential. 
In order to identify the most important contribution to this momentum dependence figure 7 shows the momentum 
dependence of the single-particle potential in nuclear matter including the various components of the M3Y interaction. 
As our discussion of the local density approximation above demonstrated that the main ingredients of the energy- 
dependence are already obtained for nuclear matter, we will concentrate our discussion on nuclear matter. From 
the analysis displayed in figure 7 it is obvious that the momentum dependence is dominated by the components of 
medium range in the exchange part of the M3Y interaction. The pion-exchange contribution is much weaker and 
the repulsive short-range components increase the attraction for larger momenta. This medium range component in 
Eq.(^ is mainly due to Vr and Var- Translated into the language of mesons this means that it is due to the exchange 
of an isovector meson with a mass of around 500 MeV. As such a meson does not exist, we must consider this as an 
effective meson to parameterize the effects of correlations or other features contained in the G-matrix. 

It is very difficult to decide whether such strong exchange effects required in the M3Y or other non-relativistic 
parameterizations are realistic or just an artefact of this parameterization. The only information available from NN 
scattering data and microscopic calculations of NN interactions is the information about the sum of direct and exchange 
term and therefore the separation into a direct and an exchange term is to some extent model-dependent. There are, 
however, some indications that such strong exchange effects do not show up within the relativistic framework. 

It is worth noting that in relativistic parameterizations of the G-matrix, the self-energy is dominated by the direct 
Hartree-contributions. This is valid for the Dirac-Hartree-Fock parameterization of po|] , which we have discussed so 
far. The same is also true for the more complete parameterization of the G-matrix derived from the Bonn potential of 
pl[ and for the independent analysis of the group in Groningen [ pT|j4^ ] . The Fock- or exchange-terms yield corrections 
to the Hartree terms, which are as large as 30 percent of the Hartree term. However, one does not meet a situation 
that the exchange terms are more important as it is the case in the non-relativistic framework discussed above. 

This is also reflected by the fact that the momentum or energy-dependence of the SEP derived from the Dirac- 
Brueckner-Hartree-Fock self-energy in nuclear matter is dominated by the energy dependence displayed directly in 
Eq . (po|) . This is reflected in figure 8, where the momentum dependence of the SEP in nuclear matter is displayed with 
(solid lines) and without (dashed lines) taking into account the momentum dependence of the Dirac self-energy. Both, 
for the Bonn potential and for the independent analysis of the Groningen group, one observes that the momentum 
dependence of the Dirac self-energy, which reflects the Fock exchange effects, is a small correction and tends to reduce 
the momentum dependence of the SEP. 



V. CONCLUSIONS 

Results for the real part of the optical potential for nucleon - nucleus scattering at low energies (E < 100 MeV) 
are discussed within a non-relativistic framework and the relativistic Dirac-Hartree-Fock approximations. The in- 
vestigations are based on relativistic and non-relativistic parameterizations of a G-matrix derived from realistic NN 
interactions. Both approaches lead to local energy-dependent potentials, which, employed in a Schroedinger equation, 
reproduce the bulk features of the energy dependence of the empirical data for nucleon scattering on and '^'^Ca. 
However, for the same density distribution of the target nucleons different local equivalent potentials are obtained. 

A local density approximation (LDA) contains the main features of the energy dependence of the optical potential. 
The LDA fails to reproduce the precise shape of the potentials leading to smaller radii. Also it does not provide 
accurate results for the binding energy. 
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Both kinds of approaches yield a similar energy- or momentum dependence. Therefore, analogous to the empirical 
studies of nucleon nucleus scattering, also the present investigation cannot give a definite answer to the question 
whether relativistic features are needed to describe these processes. The origin of the energy dependence is quite 
different in the two approaches. Within the non-relativistic scheme strong Fock exchange terms are required to obtain 
an energy dependence of the optical model as required from the empirical data. In the relativistic framework the 
energy dependence is due to the reduction of the Dirac equation to a Schroedinger equation and Fock exchange terms 
play a minor role. Therefore one may understand the strong exchange effects of the non-relativistic approach as a 
tool to incorporate the relativistic features. 

It should be useful to study the optical potential for the elastic scattering of heavy ions to investigate whether 
the empirical data for these systems allow to distinguish between the relativistic and the non-relativistic approach. 

This investigation has been supported by the "GraduiertenkoUeg Struktur und Wechselwirkung von Hadroncn und 
Kernen" (DFG, Mu 705). 
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TABLE I. Parameterization of the diflferent central terms of the M3Y potential (see Eq.(l)), using the parameterization of 
Eq.(2). Note that for the err channel the pion contribution of Eq.(3) has to be added. 



k Hk [fm-i] C<'=' [McV] C?^ [McV] Ci'=' [McV] [McV] 



1 4.0 7999.00 -4885.50 -2692.25 -421.25 

2 2.5 -2134.25 1175.50 478.75 480.00 
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FIG. 1. Density distribution and optical potential for O. The density distribution displayed in the upper part of the figure 
has been obtained from the Dirac Hartree calculation of Ref.[20]. For this distribution the local equivalent optical potential 
of nucleon scattering have been calculated for energies of 5 MeV and 100 MeV, using the Dirac Hartree (solid line) and the 
non-relativistic folding potential calculated for the M3Y potential (dashed curve) 



FIG. 2. Density distribution and optical potential for Ca. The lower part of the figure also displays the local density 
approximation to the optical potential defined in Eqs.(24) and (20). Further details see figure 1 



11 



FIG. 3. Phase shifts for / = scattering on 'O at various energies of the scattered nucleon. Results of the Dirac Hartree 
calculation are displayed in the upper part while those obtained from the non-relativistic M3Y potential are given in the lower 
part of the figure. Beside the results of the consistent calculations (solid lines) also phase-shifts obtained for the equivalent 
local potential at a fixed energy of 5 MeV (dashed line) are presented. For comparison also the phase-shifts obtained from the 
empirical fit of Ref.[3] are displayed as "experimental" data. 



FIG. 4. Phase shifts for i = scattering on Ca at various energies. Further details see figure 3. 
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FIG. 5. Calculated I — phase shifts for scattering on Ca. Results are displayed for the Dirac-Haxtree, the 
Dirac-Hartree-Fock and the local density approximation for the Dirac-Hartree approach. 



FIG. 6. Contributions to the binding energy of nuclear matter at various densities calculated from the direct part of the 
M3Y interaction (Hartree) and from the direct plus exchange part (Hartree-Fock). The kinetic energy is not included. 
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FIG. 7. Momentum dependence of the single-particle potential in nuclear matter at saturation density {kp = 1.36 fm~^) 
calculated for the M3Y potential. Results are displayed including the total exchange part of the interaction (solid line), the 
medium and short-range part (dashed line) and the short-range part only (line with circle labels). Displayed is the energy 
difference to the potential energy at momentum A; = in each case. 



FIG. 8. Momentum dependence of the Schroedinger equivalent potential (SEP) calculated according to Eq.(20) in nuclear 
matter at saturation density (kp = 1.36 fm"'^). Results are displayed which account for the momentum dependence of the 
scalar and vector terms in the self-energy (solid lines) and those for which the scalar and vector terms calculated at momentum 
k — kp were taken for all momenta (dashed lines). Results derived from the Bonn A potential (Ref. [20]) are compared to 
those of the Groningen group (Ref.[17], lines with circle labels). Displayed is the energy difference to the potential energy at 
momentum A; = in each case. 
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